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We address the computation of the determinant for the Dirac operator corresponding to a quark propagating in a back- 
ground gluon field of the following type: the gauge field is covariantly constant and self-dual inside a hypersphere and 
with quark fields satisfying bag-like (chiral- violating) boundary conditions at the surface of the hypersphere. We find that 
the parity odd part of the logarithm of the determinant corresponds to the Abelian anomaly. However the parity even 
part also depends on the chiral angle associated with the boundary. This unexpected dependence is discussed. 



1 Introduction 

The profound physical phenomena expected to arise 
from QCD, such as confinement, spontaneous chiral 
symmetry breaking and consequences of the Abelian 
anomaly, have inspired a range of simplified mathe- 
matical settings for the QCD equations and searches 
for exact solutions to such problems. The exact self- 
dual and finite action solutions to the Yang-Mills clas- 
sical equations known as the instanton are one exam- 
ple of this. In the quark sector, a similar endeavour has 
developed involving searches for exact solutions to the 
Dirac eigenvalue spectrum for given nontrivial back- 
ground gluon fields, and in turn computations of the 
quark determinant (and thus the quark free energy) in 
such gluonic backgrounds. Often using zeta function 
regularisation or heat-kernel methods, such problems 
have now been solved for various compact spaces in 
various dimensions with and without nontrivial gluon 
fields in the background. In this paper we shall address 
such a computation for the quark determinant in what 
we have previously called a "domain-like" gluon back- 
ground field, based on our earlier solution to the Dirac 
eigenvalue problem for such a configuration. 

The elements of the problem are as follows: we con- 
sider a four-dimensional Euclidean hypersphere of ra- 
dius R. On the surface, fermion fields satisfy a chi- 
rality violating boundary condition reminiscent of bag 
models of the nucleon, but now for four-dimensional 
domains with nontrivial field internal to the domain. 
In various space-time dimensionalities, recent works 
have addressed such boundary conditions [ 121 EH Bl 
|S1 E] illustrating the relevance of the problem to the 
realisation of chiral symmetries in QCD. For exam- 
ple, in the fermion boundary condition is regarded 
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as a mechanism for dynamical generation of the CP- 
breaking 0-term. 

An interesting feature of these boundary conditions is 
that they render the spectrum of the Dirac operator 
asymmetric under A — > — A which in turn has fascinat- 
ing consequences. In particular, a careful analysis of 
the fermion determinant leads to an additional asym- 
metry spectral function [|HHHIE3]> as well as t ne usual 
zeta function type contributions. In this respect, the 
four dimensional theory we tackle resembles three di- 
mensional field theories where parity violating effects 
are well-known and accessed through such asymme- 
tries [ 8 . The new element we introduce here is the 
presence inside the hyperspherical region of a non- 
vanishing covariantly constant, self-dual gluon field. 

This type of gluon background field arises in the "do- 
main model" ansatz for the QCD vacuum [^3> con ~ 
sisting of an ensemble of such hyperspheres with ran- 
domly assigned field orientations and self-duality or 
anti-self-duality. The fermion boundary condition in 
this case arises as a result of requiring finite action 
fluctuations in the presence of a singular field defining 
the boundary. A chiral angle associated with fermion 
boundary condition here becomes an additional ran- 
dom variable associated with a given domain, such 
that averaging over the angle for an ensemble of do- 
mains restores chiral symmetry for the ensemble [I12|. 
In such a model, we have previously demonstrated 
confinement of static quarks anc l have shown that 
the chirality of low lying Dirac modes exhibit strong 
correlation with the duality of the gluon field inside the 
domain [^j reminiscent of lattice QCD calculations [ 
1131 1141 IT5j and consistent with the dominance of the 
chiral properties of the lowest modes in condensate for- 
mation, as expressed in the well-known Banks-Casher 
relation [llbj. The computation of the quark determi- 
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nant, tackled in this paper, is crucial for studying chi- 
ral symmetry realisation in the domain model, and can 
be important for understanding this problem in QCD 
itself. In particular, the determinant is precisely the 
relevant quantity for obtaining the Abelian anomaly 
in the path- integral approach [ I17j . The calculation 
presented here recovers the Abelian anomaly, which is 
identified as the parity odd part of the logarithm of de- 
terminant. The anomaly gets contributions from both 
the standard ^-function of the squared Dirac opera- 
tor and the spectral asymmetry function, since in the 
presence of chirality violating boundary conditions the 
spectral ^-function and spectral asymmetry function 
do not have direct physical meaning as parity even and 
odd respectively. The dependence of the anomaly on 
the chiral angle turns out to be multivalued, which is 
to be expected since the topological charge here is not 
restricted to integer values. The parity even part also 
shows a dependence on the chiral angle (the parameter 
of boundary conditions), which we consider to be more 
an artifact of the incompletness of our calculation than 
an established property of the determinant. Here we 
will only identify those potential contributions which 
are beyond the scope of the present calculation, and 
which can eliminate the chiral angle dependence of the 
parity even part. 

2 Boundary Conditions, Dirac opera- 
tor and eigenvalues 

Our starting point is the eigenvalue problem for the 
Dirac operator subject to the bag-like boundary con- 
dition 

pip(x) = Xip(x), (1) 
fj(x)j 5 e^ip(x) = i/>(x), x 2 = R 2 . (2) 

Here rj^(x) = x^/\x\, is the covariant derivative in 
the fundamental representation, 

-Dp = dp- iBf, 
h = t a n a , 

where the (anti-)self-dual tensor is constant, and 
the Euclidean 7— matrices are in an anti-hermitean 
representation. To be specific, we will take the field 
as being self-dual. The magnitude of the topological 
charge per domain will be a useful quantity for later, 
and here it turns out to be [|n] just q — B 2 R 4 /16. 
As in [[I], eigenvalues defined by Eqs. JIJ and J2J) are 



real for all real *&. In this case the boundary condition 
for ip required by the condition 

i>{x) r}(x)ip(x) = 0, x 2 = R 2 , 

is given by hermitian conjugating Eq. ft}. Complex 
values of the chiral angle require in general construc- 
tion of a bi-orthogonal basis. Crucial to the solution 
derived in [ I12| was the decomposition 

i> = iffX + <P> $ = iX fi + <P, (3) 

ip and x having the same chirality. Solutions are then 
labelled via the Casimirs and eigenvalues 

K? = Kl^^ + 1), fc = 0,l,...,oo 

m h2 = -k/2, -jfc/2 + 1, . . . , fc/2 - 1, fc/2, 
k = 0,1,2,..., 

corresponding to the angular momentum operators 
K1.2 = |(L ± M) with L the usual three-dimensional 
angular momentum operator and M the Euclidean 
version of the boost operator. Thus the solutions for 
the self-dual background field can be written 

^fcm 1 =*^ l +^ 1 . (4) 

The minus sign label on these solutions refers to their 
negative chirality, 75^ = —<p, 75% = — Xj an d is the 
only choice of chirality for which the boundary con- 
dition Eq. J2J can be implemented for the self-dual 
background field [ 12: . (For an anti-self-dual field, the 
chiralities become positive.) The label k = ± refers to 
the polarisation of the spinors according to the com- 
bined colour-spin projection specified by the projector 

K = 7V+£ K + A^£_ K (5) 

with 

N ± - i(l±n/|n|), 
E± = i(l±E.B/S) 

being respectively separate projectors for colour and 
spin polarizations. 

In the ^-regularised determinant an arbitrary scale fi 
appears and it is convenient to work with scaled vari- 
ables 

(3 = \Z2§/fx, p = nR, £ = X(B)/fi, & = A(0)/m, (6) 
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and where the dimensionless quantity z — BR 2 /2 = 
f3 2 p 2 /4 appears prominently. Then, as derived in [ I12j . 
the equations determining the (rescaled) eigenvalues £ 
for the two polarisations that appear in the problem 
are 



M(-e/(3 2 ,k + l,z) 

■d sP A/r /i , :„> ,2 



+e 

+ 2(k + l] 
for k = — and 



M{l-^/l3\k + 2 lZ ) = 



(7) 



function respectively: 



<q(s) 



k,7i,< 



,2s 



[AL(0)] 2s 



, ^ / sgn(A L (5)) 
Bgn(Ag„(0)) 

|AL(o)l s 



(11) 



(12) 



M(fc + 2-£ 2 //3 2 ,fc + 2,z) 
+e ^[M(fc + 2-e 2 //3 2 ^ + 2,z) 



fc + 2-£ 2 //3 2 
fc + 2 



M(/c + 3-C 2 //3 2 ; fc + 3,z) 



(8) 



= 



for k = +. The confluent hypergeometric, or Kummer, 
functions M(a, b, z) emerged originally as solutions of 
the radial equation for a scalar field in the presence 
of the above covariant constant gauge field on a finite 
hypersphere [ITT], 



3 Computing the determinant 

We now consider the quark determinant in a single 
self-dual domain of volume v = ir 2 R 4 /2 



detf 



ip 



exp 



fc,n,mi 

exp{-C'(s)} s=0 



0) 



where the arbitrary scale p, plays a role similar to the 
renormalization point in other regularisation prescrip- 
tions. The zcta function £(s) breaks up into two parts [ 
|S] respectively symmetric (S) and antisymmetric (AS) 
with respect to A — ► —A 



C(*) = Cs(s) + Cas(s), 



(10) 



It should be stressed that in the presence of bag-like 
boundary conditions £g an d Cas do not have the mean- 
ing of parity conserving and parity violating terms 
since a parity transformation in terms of eigenvalues is 
given by A(i9) — » — A(— i?) and both spectral functions 
contain parity conserving and violating terms. Thus 
the determinant for a given parameter i9 is defined by 



C / (O) = Q0 2 (O)±i^(O)Tifr ? (O)V 



(13) 



with the normalization chosen in Eq.Q such that 
C'(0) vanishes as B — > 0. 

Spectral sums can be computed using a representation 
of the sum as a contour integral of the logarithmic 
derivative of the function whose zeroes determine the 
spectrum, 



y- 



1 

27Ti 



(14) 



see for example [EIL where the zeroes of /(£) = are 
£ = A and the contour is chosen such that all zeroes are 
enclosed. With real parameter the poles lie on the 
real axis, and there is no pole at the origin for any 
By deforming the contour and accounting for the van- 
ishing of contributions to the integral at infinity, the 
expressions arising from Eq. (|14f> can be transformed 
into real integrals. The following representations for 
the two spectral functions are eventually obtained 



with 



Cas(s) = ^{l-e^ s ) V (s) 



and the key quantities devolve into the zeta function 
for the squared Dirac operator and the asymmetry 



,(*) 



rj(s) 



, sin(7rs) 



E fcl ~ 



fe=i 

x / £|*<mm), 



, cos(7rs/2) 



k=l 



/ -$(fc,tz,i?), 

Jo t s dt 



(15) 



(16) 
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with 'J' and <f> being the sum of contributions from the 
two polarisations, taking the form 



Riemann zeta function. The resulting £^>2 (s) has the 
structure 



$(k,t\z,i}) 

<$>(k,t\z,&) 

and where 

A-(k,t\z) 
B^{k,t\z) 

A+(k,t\z) 
B+(k,t\z) 



K=± 

E 1 * 



K=± 



A 2 K (k,t\z) + e 2 «BZ(k,t\z) 
A 2 {k,t)+e 2i} B 2 {k,t) 

A K {k,t\z) + te i) B K (k,t\z) 
A K (k,t\z) - ie^B K (k,t\z) 



M(—r^-,k + l,z 



4z 



2+2 „2 



ktp „,, k 2 t 2 p 2 , 
' M(l + —^-,k + 2,z), 



2(fc + l 



43 



M( _^,fc + l,-z), 



2z 
ktp 

k + 1 



4z 

/r 2 f 2 n 2 

M(-^4^-,Jfe + l,-z) 



4z 



xM( 



fe + 1 

fc 2 tV 



A(fc,t) 
B(k,t) 



4z 



,fc + 2, -z) 



(17) 



The Bessel functions Ik emerge from the limit B — > 
with the normalization of the determinant as specified 
above. Analytical continuation of (p2 (s) and rj(s) can 
be done by means of standard methods (for an exten- 
sive review see [El)- These functions can be related 
to functions such as the Riemann zeta function Cr( s ) 
with known analytical properties. Our particular ex- 
ample of such a calculation is clearly technically quite 
complicated. First of all, expansions of the confluent 
hypergeometric functions in 1/k, similar to the Debye 
expansion of Bessel functions, is required, 



M{—Z-,k + l,z) 



4z 



C(tp, k) 



M n (tp,z) 



n=0 



(18) 



The form of the prefactors and the M n (x, z) functions 
will be presented in detail elsewhere [E, • It suffices to 
say that the prefactor C(x, k) cancels out of $ and W 
after which the procedure is straightforward if some- 
what tedious: the expansions Eq. (|18fl are inserted 
in £(s) and r)(s) so that the integrals over t can be 
evaluated term by term. After exchange of order of 
summation over k and n (which is the most subtle 
step here), sums over k can be read off in terms of the 



, 2 ( s ) = sp 2 S £mM g +n _ 2 



n=0 



+6{f>i(s), (19) 



where the term in the second line denotes those contri- 
butions coming from interchange of the order of sum- 
mations over n and fc, which are potentially present. 
A similar structure appears for the asymmetry spec- 
tral function as well. Finally we are interested in 
the decomposition of this expression around s = 0. 
The result for the first line in Ea. (|19|l comes from the 
n = 2 term alone and can be calculated with relative 
ease since only the lowest coefficients Mi and Mi in 
Eq. (|18f) contribute. However, as occurs in even sim- 
ple problems the second term in Ea. (|19l) is much 
more difficult to compute. To achieve this one needs to 
know coefficients f(z 2 \n) as a function of continuous 
variable n. Below we will concentrate on the contri- 
bution given in the first line alone, bearing in mind 
the necessity of a complete analysis. The final results 
for the first term in Eq. Ijl9|l and its analogue in the 
asymmetry function ?y(s) are then summarised in the 
following equations: 



O(o) 


= o, 






C>(0) 


= z 2 


-i-ln2 + ln(l + e 2 ") 




77(0) 


z 2 
2 


z 2 

H Arctan(sinh(tf)) 

7T 


(20) 



and hence 



,2 r 



C'(o) = - 



1 IT 

In2±i- 

4 2 



(21) 



+ In (1 + e™) T iArctan(sinh(tf))] + <5C'(0), 

where <5C'(0) denotes contributions potentially coming 
from the above-discussed non-commutativity of sum- 
mations, which have been not calculated here. For 
an anti-self-dual background field, we merely have to 
perform a parity transformation on the fields. For the 
spectrum it involves A($) — > — A(— For the deter- 
minant, one need only take $ — > — $ and rj(s) — > —r](s) 
independently. 

4 Discussion and Conclusions 

For interpretation of the result Eq. I|21(l one needs to 
continue $ to imaginary values. The appropriate con- 
tinuation is i? — > ia + iir/2. In this case, one can 
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cleanly separate parity odd and even parts since for 
the self-dual field one gets 



C(0) - y 



-- + ln(l ± cos(a)) + i(a + ttv) 



with v G Z. while the anti-self-dual field leads to 



cm - y 



-- + ln(l ± cos(a)) — i(a + -kv) 



The imaginary part is given in a form which mani- 
fests the 27T periodicity in a. Recalling that the ab- 
solute value of the topological charge associated with 
the gluon field in a domain is q = z 2 /4, the parity odd 
part becomes 

± 2iq{a mod if) 

with the sign correlated with the duality of the field, 
and can be regarded as the axial anomaly. This agrees 
with the estimation of obtained in a more general 
approach, up to a term ±iqir. As far as we can see 
this difference occurs due to contributions in our cal- 
culation from both £™ 2 (0) and the asymmetry spec- 
tral function 77(0), which provides an extra ±iqir to 
give the above anomaly term. The contribution of the 
asymmetry spectral function was not considered in [ 

Final conclusions about the validity of the chiral angle 
dependence of the parity even part, as well as an anal- 
ysis of the strong background field limit, necessarily 
require computation of the potentially nonzero term 
6( p in Eq.jfnj- 



3. M. Bordag, E. Elizalde, K. Kirsten, J. Math. 
Phys. 37 (1996) 895. 

4. J.S. Dowker, J.S. Apps, K. Kirsten, M. Bordag, 
Class. Quant. Grav. 13 (1996) 2911. 

5. C.G. Beneventano, M. De Francia, 

E.M. Santangelo, Int. J. Mod. Phys. A 14 (1999) 
4749. 

6. K. Kirsten, Nucl. Phys. Proc. Suppl. 104 (2002) 
119. 

7. V.N. Marachevsky, D.V. Vassilevich, 
|hep-th/ 0309019 (2003). 

8. S. Deser, L. Griguolo, D. Seminara, Phys. Rev. 
Lett. 79 (1998) 1976. 

9. E. Elizalde, M. Bordag, K. Kirsten, J. Phys. A 31 
(1998) 1743. 

10. C.G. Beneventano, E.M. Santangelo, A. Wipf, 
hep-th/0205199 (2002). 

11. A.C. Kalloniatis, S. N. Nedelko, Phys. Rev. D 64 

(2001) 114025. 

12. A.C. Kalloniatis, S.N. Nedelko, Phys. Rev. D 66 

(2002) 074020. 

13. I. Horvath, N. Isgur, J. McCune, H.B. Thacker, 
Phys. Rev. D 65 (2002) 014502. 

14. R.G. Edwards, U.M. Heller, Phys. Rev. D 65 
(2002) 014505; T. Degrand, A. Hasenfratz, Phys. 
Rev. D bf 65 (2002) 014503; I. Hip et al, Phys. 
Rev. D 65 (2002) 014506; T. Blum et al, Phys. 
Rev. D bf 65 (2002) 014504. 

15. C. Gattringer et al, Nucl. Phys. B 618 (2001) 
205; C. Gattringer et al, Nucl. Phys. B 617 
(2001) 101. 

16. T. Banks, A. Casher, Nucl. Phys. B 167 (1980) 
215. 

17. K. Fujikawa, Phys. Rev. D 21 (1980) 2848; 
Erratum-ibid. D 22 (1980) 1499. 

18. E. Elizalde, Ten Physical Applications of 
Spectral Zeta Functions, Lecture Notes in Physics, 
m35, Springer- Verlag (Berlin, Heidelberg, New 
York) 1995. 

19. A.C. Kalloniatis, S.N. Nedelko, in preparation. 



Acknowledgments 



S.N.N, was supported by the DFG under contract 
SM70/1-1 and, partially, by BFBR grant 01-02- 
1720. We acknowledge the hospitality of the Institute 
for Theoretical Physics III, University of Erlangen- 
Nuremberg where substantial parts of this work was 
done. A.C.K. is supported by a grant of the Australian 
Research Council. 

References 

1. A. Wipf and S. Duerr, Nucl. Phys. B 443 (1995) 
201. 

2. S. Duerr, A. Wipf, Ann. Phys. 255 (1997) 333. 



